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We report on the effects of quantum interference induced by transmission of an arbitrary number 
of optical quantum states through a multiple scattering medium. We identify the role of quantum 
interference on the photon correlations and the degree of continuous variable entanglement between 
two output modes. It is shown that the effect of quantum interference survives averaging over all 
ensembles of disorder and manifests itself as increased photon correlations giving rise to photon 
anti-bunching. Finally, the existence of continuous variable entanglement correlations in a volume 
speckle pattern is predicted. Our results suggest that multiple scattering provides a promising way 
of coherently interfering many independent quantum states of light of potential use in quantum 
information processing. 

PACS numbers: 42.25.Dd, 42.50.Lc, 78.67.-n 



Studies of wave propagation in disordered media have 
revealed a range of fascinating wave phenomena, includ- 
ing Anderson localization {l|, enhanced coherent back 
scattering {3| , and universal conductance fluctuations Q . 
These phenomena, originating from wave interference, 
appear after averaging over all configurations of disor- 
der, in the mesocopic regime [1,0]. Mesoscopic effects are 
common to a large class of elastic wave problems, includ- 
ing transport of electrons through artificial nanostruc- 
tures, advancement of sound waves in turbulent fluids, 
and light propagation through complex dielectric media. 

While much attention has been devoted to propaga- 
tion of classical light waves through random media over 
the years, the influence of the quantum nature of light 
has only recently become an active field of research. The 
development of a theoretical framework to handle mul- 
tiple scattering in quantum optics 0] triggered the 
interest in understanding quantum optical properties of 
disordered media [6-16]. Studies include quantum noise 
properties d, absorbing or amplifying media [tI, 
degradation of polarization entanglement [lo| . and spa- 
tial photon correlations ^12nl6i l. Recently it was shown 
experimentally that light-matter interaction is strongly 
enhanced in disordered photonic crystal waveguides, en- 
abling cavity quantum electrodynamics with Anderson- 
localized modes [l7|- Optical quantum information pro- 
cessing schemes rely on interference among multiple inde- 
pendent quantum states, i.e. quantum interference (QI), 
to generate quantum correlations and entanglement. The 
possibility of using multiple-scattering media to interfere 
independent quantum states is appealing since it is inher- 
ently scalable to multiple input states. To this end, meso- 
scopic quantum interference effects that persist even af- 
ter averaging over all ensembles of disorder, would be re- 
quired in order to obtain robust and predictable quantum 
correlations. So far, interference of quantum light in a 
scattering medium has only been described in ID random 
walk models [l5j | or for diffusive transport where interfer- 
ence effects wash out after ensemble averaging 12, lif . 



FIG. 1: (color online). Sketch of propagation through a dis- 
ordered waveguide of length L and transport mean free path 
£. Quantized light is incident to the left and the correlations 
between two output modes on the right are analyzed. The 
operators ai and Sq corresponds to the creation operators of 
modes i and a, where Roman and Greek subscripts denote in- 
put and output modes respectively. The correlations between 
the two different output modes a and /3 are analyzed. 



In this Letter, we investigate QI induced by combining 
an arbitrary number of independent quantum states in 
a random multiple scattering medium in the mesoscopic 
regime. We identify the role of QI on the degree of pho- 
ton number correlations between two transmission paths 
through the medium and the degree of continuous vari- 
able entanglement. Surprisingly QI of photons is found 
to survive after averaging over all configurations of disor- 
der in the mesoscopic regime, i.e. the induced quantum 
correlations have deterministic character despite the un- 
derlying random multiple scattering processes. At last 
we discuss the feasibility of experimentally verifying our 
theoretical predictions. 

Let us now introduce the model for propagation of 
quantized light through a linear, elastic, multiple scatter- 
ing medium of length L and transport mean free path i, 
see Fig. [TJ We apply the scattering matrix for the prop- 
agation of light and use random matrix theory on the 
scattering elements. The approach describes effectively a 
quasi-lD model of an A'^-mode waveguide, but is known 
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FIG. 2: A diagrammatic representation of three terms of the 
form t'^it*pjti3ktai{alajakai). Diagrams (a) and (b) involve 
only intensity transmission of the input modes while (c) shows 
quantum interference between two input states. The three 
diagrams are the only ones that survive ensemble averaging. 



also to accurately predict propagation in 3D slab geome- 
tries Q. We relate the photon annihilation operators 
CLa (oi) of output (input) modes a (i) by cLa=J2itaiai, 
where the summation is over all N possible input modes 
at each end of the waveguide and tai denotes the complex 
scattering matrix element. Experimentally such a system 
could, e.g., be realized in titania powder samples |14j or 

disordered photonic crystal waveguides [13, m. 

As a measure of QI, we introduce the 2-channel photon 
correlation function 



Cap 



where 



Afiafiis = {fiafin) - {na){np). 



(1) 



(2) 



The brackets denote quantum mechanical expectation 
values and na=a^aaa is the output photon number oper- 
ator. The degree of entanglement is quantified in terms 
of the quadrature variance product (QVP) 



addition of the intensities associated with the different 
propagation paths through the medium, as it is the case 
for the diagrams (a) and (b). The latter gives rise to 
QI between the input states. Diagram (c) is an example 
that show such interference. We note that all additional 
diagrams not shown in Fig. [2] are interference diagrams. 

For diffusive transport the intensity transmission coef- 
ficients are exponentially distributed while the phase is 
uniformly distributed [20[. Consequently, this allows us 
to generate a set of normalized amplitude transmission 
coefficients valid for a single realization of disorder. To be 
specific, we choose a waveguide of A^^IO^ so that the av- 
erage single channel transmission is T=g/N^=l/300 [U], 
where g is the normalized average conductance. The 
modes are written in a 10 by 10 grid to represent spa- 
tial wavevectors in the transverse plane {kx, ky). We 
evaluate the 2-channel photon correlation function us- 
ing Fock states \n) as input. Illuminating only a sin- 
gle input channel with a two-photon Fock state |2), 
Cq,^=— 1/2 for all modes a and /3, independent of the re- 
alization of disorder. For single photons incident in two 
different input modes |1,1) the spatial photon correla- 
tions fluctuate between -1 and 0, see Fig. [3] (a). This 
is a manifestation of QI in a speckle pattern as ob- 
served in Ref. We furthermore evaluate Sa/s for 
two quadrature-squeezed input states \(^i)=exp[^Q — 
^Ci(fl^)^] |0) where Ci=\Ci\e-'''^' contains squeezing ampli- 
tude. Id, and phases, 0]. We investigate two or- 
thogonally oriented squeezed beams, i.e. (t>i=0 and (t>j=Tr 
and chose |C|=0.15 corresponding to experimentally ob- 
tainable parameters [3]. Fig. [3] (b) is a calculation of 
Eap for different output modes and show that entangle- 
ment (£c[/3<l) can be induced by multiple scattering. If 



(3) 



where Xa—^{d.l,+a.a) and Ya=^{al^—aa) are quadra- 
ture operators. The QVP determines the ability to pre- 
dict a measurement in mode f3 given the result of a 
measurement on mode a, and for ea/3<l (>1) the out- 
come is predictable below (above) the quantum noise 
limit. Here, eap<l implies that the quantum state of the 
two output modes a and /? is unseparable, i.e. entan- 
gled [19|. Eqs. dl]) to ([3]) can conveniently be evaluated 
diagrammatically, by representing the propagator taicii 
by an arrow connecting input mode i to output mode 
a. Since the scattering matrix is unitary, i* joj repre- 
sents the time-reversed path. E.g. considering two input 
modes, evaluating Eq. ([2]) yields 2* different terms of the 
form t'^it*^jtj3ktai{(ii(ijO,kO,i)- Three typical diagrams are 
shown in Fig. [5J As an example, the contribution from 



diagram (a) is \tai\^\tpi 



:), where 



denotes 



normal ordering and ni=a\ai is the input photon number 
operator. The diagrams can be classified into intensity 
and interference diagrams. The former is an incoherent 
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FIG. 3: (color online), (a) The 2-photon correlation Ca/s for 
two single-photon Fock input states, |1, 1), showing large fluc- 
tuations due to quantum interference, (b) The degree of en- 
tanglement logj^Q(ea/3) for the two quadrature squeezed input 
states described in the text, where the gray-scale corresponds 
to non-entangled states (logjQ(ea/3)>0) while the colored ar- 
eas with crosses display the entangled states (logjQ(ea/3)<0). 
(a) and (b) are obtained with the same realization of dis- 
order for diffusive transp ort, w here t he ph ase is random 
and P(|fciiP)=exp( — Itaip/ltaiP) with \tai\^=r the average 
single-channel transmission. 
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we change the squeezing phases of the input states the 
modes that display entanglement change. From knowl- 
edge of the transmission matrix one could thus specify 
the mode in which the entanglement should occur by 
changing the phases of the squeezing parameters. This 
could potentially be achieved with the recent scheme to 
measure the complex transmission matrix for light prop- 
agation through a disordered medium [23| . The inher- 
ent ability of a multiple scattering medium to mix many 
modes shows the scalability of the approach of potential 
use in quantum information processing. 

Let us next consider the effects of QI after ensemble av- 
eraging. For this we need the averaged amplitude trans- 
mission coefficients which are given by [23 | 



f* f 



(4a) 

iai^pjtpktai = {Ci5a5jk + C2Sik6ji) , (4b) 

with Ci and C2 the short and long-range correlation func- 
tions respectively and the bar denoting ensemble averag- 
ing. When carrying out the ensemble average only loop- 
type diagrams survive, i.e. only those shown in Fig. [2] 
The values of diagrams (a) and (b) are proportional to 
C1+C2 and Ci respectively while diagram (c) is propor- 
tional to C2. This can be intuitively understood since 
Ci expresses the intensity fluctuations in a single speckle 
spot, while C2 is related to the intensity correlations be- 
tween two speckles ^] . Diagram (c) displays the QI term 
that survives ensemble averaging. From the values of 
C2 and the normalized average conductance g, we define 
the transitions from the quasi-ballistic to the weakly dis- 
ordered regime (6*2=0) and from the weakly disordered 
to the localized regime (g— 1). The mesoscopic regime 
is defined as the regime in which two speckle spots are 
correlated after ensemble averaging (C2>0). The values 
of Ci and C2 depend on the number of modes N and 
the decree of disorder, which is contained in s=L/£ and 

We define the ensemble-averaged 2-channel correlation 
function as 



C. 



a/3 = 



{na){nfj) 
{C1+C2) 



\ i / i 



{n^)) 



(5) 



-1. 



C2\J2{n.r + 2j2\{4&s)\' 



i,j>i 



The ensemble averaged QVP is 
£a/3 =1 + 4t^ AajcLi 



(6) 



■At' 



Ci i Aaja, j + C2^ AajcijAajc 




Normalized Lengtli s = j 



FIG. 4: (color online). The ensemble averaged 2-channel 
photon correlations, C^/? and the degree of entanglement 
logjg(eQa) versus s—L/£ for A'^ = 50. Solid blue curves show 
for Fock input states with a total of two photons and 
dashed red curves are for three photons. The difference in Cap 
between having one and more input states is due to quantum 
interference. The green dash-dotted line shows logj^Q(eci;9) for 
the two quadrature squeezed input states described in the 
text. The vertical line at s—2 indicates the crossover from the 
quasi-baUistic to the mesoscopic regime (C2=0) and the one 
at s«40 shows the transition to the localized regime (<?=!). 
The symbols o, x, and * on the abscissa correspond to the 
experimental structures studied in Refs. [l^, [ij, and [T^. 
as discussed further in the main text. 



In Fig. |4]Cq^ and logig(ea^) are plotted versus s. For 
Fock input states the s dependence of Cq,/? is a direct 
measure of QI since disregarding the QI terms implies 
that Cap only depends on the total number of input pho- 
tons. First consider having two photons incident in only 
one mode, |2), Ca/3=— ^ independent of s. With the pho- 
tons in two different input modes, |1, 1), the correlations 
on the contrary depend on s. Only at the transition to 
the mesoscopic regime we have Ca^=— ^. This value cor- 
responds to the correlation between two equally probable 
output modes for two classical non-interacting particles. 
This is due to two simple reasons: (i) the ensemble aver- 
aging makes all output modes equally possible and (ii) in 
this limit transport is diffusive and thus all interference 
effects are washed out. As disorder is increased, Cq./? in- 
creases signifying that the probability that the two pho- 
tons arrive at two different positions increases although 
remaining anti-correlated (Cct/3<0). The increased cor- 
relations saturate in the localized regime since C2 tends 
towards Ci [2^ . The variations in Ca/s can be attributed 
to QI amongst the input channels, which causes the pho- 
tons to antibunch. For three photon Fock states the more 
incident modes the larger the effect of QI on Ca^- For 

1,1,1), Cap tend towards zero in the localized regime 
signifying that the output modes become uncorrelated. 
If the number of input modes is increased further the 
output modes will become correlated (Cq^>0) in the lo- 
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calized regime. This means that detection of a photon 
in one mode on average increases the probabihty of de- 
tection of a photon in another mode, which is in striking 
contrast to the behavior of diffusive transport. Having 
single-photon states in n input modes and letting n go to 
infinity makes Cap approach unity far into the localized 
regime, which is the value obtained for thermal light. The 
variation of Sa/i with s is plotted as the green curve in 
Fig. m for the same quadrature-squeezed inputs as for the 
single realization of Fig.[3l In the mesoscopic regime both 
Ci and C2 are positive and thus eQ^>l. The value of 
approaches unity in the localized regime since the trans- 
mission decreases so that contributions from vacuum fluc- 
tuations dominate. Continuous variable entanglement in 
the transmission is therefore predicted to vanish after 
ensemble averaging. We anticipate that this might be 
different in the reflection due to coherent effects as en- 
hanced backscattering, but this is outside the scope of 
the present work. 

Finally, we address the experimental feasibility of the 
proposal. In Fig. 2] we indicate the position of three 
existing multiple scattering structures from the litera- 
ture, where the number of modes N has been scaled to 
match the value used in the calculations. Ref. [l^ con- 
cerns transmission through two scattering surfaces, which 
mimic a multiple scattering medium with s = 2. This 
corresponds to the diffusive limit where QI will be present 
in the speckle pattern but not survive ensemble averag- 
ing. In Ref. [IJ] a titania powder is used with sample 
length L=20 fiui and transport mean free path £«0.9 /xm, 
which corresponds to the mesoscopic regime of s>2. Such 
sample support a large number of modes (iV>10^) and 
thus (73>1, which means that this type of sample is in the 
weakly disordered regime where QI effects are modest, cf. 
Fig. 131 This illustrates the importance of using multiple 
scattering samples supporting only few modes in order 
to observe QI. A disordered multimode photonic crystal 
waveguide is exactly such a system and for Nfnb together 
with the typical experimental parameters of £«20 //m and 
L=100 fiui [3| gives rise to sizeable QI effects that will 
be observable in an experiment, cf. Fig. |4l 

In conclusion we have theoretically predicted that sur- 
prisingly QI of light survives multiple scattering even af- 
ter ensemble averaging and can be employed to mix sev- 
eral input quantum states of light. By calculating the 
2-channel photon correlation function in the case of Fock 
input states, we showed that the difference between one 
and several input modes for a fixed number of total pho- 
tons is a direct measure of QI. The effect was found to 
give rise to spatial photon anti-bunching in the meso- 
scopic regime and to increase with the degree of disorder. 
The experimental feasibility of the proposal was investi- 
gated based on existing multiple scattering samples from 
the literature, and we found that multimode disordered 
photonic crystal waveguides are promising candidates for 
an experimental demonstration. We furthermore investi- 



gated continuous variable entanglement induced by mul- 
tiple scattering of squeezed light that was predicted to 
show up in single realizations of disorder but to vanish af- 
ter ensemble averaging. Our work may provide a promis- 
ing new route to coherently combine many independent 
quantum states of light that is inherently scalable since 
the multiple-scattering process mixes all input states. 
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